Dynamic transition to spontaneous scalarization in boson stars 
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We show that the phenomenon of spontaneous scalarization predicted in neutron stars within 
the framework of scalar-tensor tensor theories of gravity, also takes place in boson stars without 
including a self-interaction term for the boson field (other than the mass term), contrary to what 
was claimed before. The analysis is performed in the physical (Jordan) frame and is based on a 3+1 
decomposition of spacetime assuming spherical symmetry. 
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I. INTRODUCTION 



Scalar-tensor theories of gravity (STT) are alternative 
metric theories of gravitation where a spin-0 gravitational 
degree of freedom <fi accompanies the usual tensor spin-2 
modes (see Ref. [l[ for a review). In the so-called Jordan 
frame, the scalar field <f> couples non-minimally to the 
curvature while in the Einstein frame, it couples non- 
minimally to the ordinary matter fields. Scalar tensor 
theories are perhaps the most analyzed and competitive 
theories of gravitation after general relativity, the most 
prominent example being the well-known Brans-Dicke 
theory [2[. Nevertheless, is only recently that many issues 
and theoretical discoveries concerning STT have been 
settled. In the cosmological context, STT have been pro- 
posed as alternatives to dark energy in order to explain 
the accelerated expansion of the Universe while avoid- 
ing the so called coincidence problem which is associated 
with a cosmological constant In the astrophysical 

context, Damour and Esposito-Farese 0, [l(J discovered 
that neutron star models (polytropes) within STT may 
undergo a phase transition that consists in the appear- 
ance of a non-trivial configuration of the scalar field <f> 
in the absence of sources and with vanishing asymptotic 
value. Such configurations are endowed with a new global 
quantity termed scalar charge. Due to the similarities of 
this scalarization process with the spontaneous magne- 
tization in ferromagnets at low temperatures, these au- 
thors coined the term spontaneous scalarization (SC) to 
describe this phenomenon. Further analysis |ll]4l3j con- 
firmed that this phenomenon also takes place in realistic 
neutron star models and occurs independently of the de- 
tails of the equation of state (EOS) used to describe the 
nuclear matter. 

The stability analysis for the transition to SC was first 
performed by Harada [3, HE} . It is now understood 
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that SC arises under certain conditions where the ap- 
pearance of the scalar field gives rise to a configuration 
that minimizes the star's energy (the ADM mass) with 
fixed baryon number. This interpretation can in fact ac- 
commodate the Newtonian expectations despite the fact 
that the effective gravitational constant decreases during 
the transition [l3[ . Following the analogy with ferromag- 
netism, the ADM mass plays the same role as the (free) 
energy of the ferromagnet, while the baryon mass is the 
analogue of the inverse of the temperature. The order 
parameter is the scalar charge which mirrors the mag- 
netization. One important aspect of this phenomenon 
is that it occurs even when the parameters of the the- 
ory satisfy the stringent bounds put by the solar system 
experiments. The important point is that SC appears 
precisely when the asymptotic (cosmological) value 0o 
of the scalar field vanishes. Therefore, the phenomenon 
arises even when the associated effective Brans-Dicke pa- 
rameter of the theory is arbitrary large (see Sec. |V]). On 
the other hand, the binary pulsar does put limits on the 
magnitude of SC. In some classes of STT, these bounds 
restrict the non-minimal coupling to the curvature. How- 
ever, the bounds are no so stringent since the couplings 
can still be of order unity [Io| ■ 

More recently, studies of neutron star oscillations 
within STT reveal that, in addition to the emission 
of scalar gravitational waves, SC can also disturb the 
quadrupolar gravitational radiation as compared to the 
corresponding signals in general relativity (GR). There- 
fore, even if the detection channels of scalar gravitational 
waves are "switched off" , the detection of gravitational 
waves of spin-2 coming from these sources might validate 
STT or put even more stringent bounds on their param- 
eters. Of course, the direct detection (or lack) of scalar 
gravitational waves would also help to discriminate be- 
tween several alternative theories. In this regard it is 
important to emphasize another striking feature of STT. 
While GR predicts only quadrupolar gravitational radi- 
ation in the "far zone" , STT predicts monopolar gravi- 
tational waves [16], so that even in spherical symmetry 
scalar waves can be emitted. This is in part because 
Birkhoff 's theorem does not apply in this case. The new 
polarization scalar mode is of breathing type since it af- 
fects all the directions isotropically [16| . In particular, 
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during the SC process in spherical neutron stars this kind 
of radiation might be emitted. In fact, usin g a fully rela- 
tivistic spherically symmetric code, Novak [l2j not only 
confirmed the dynamical transition to the scalarization 
state but also the emission of such scalar waves. 

In a more recent analysis, Whinnett [13] corroborated 
that SC can also occur in boson stars (see Ref. [l8| for a 
thorough introduction to the subject), but only if these 
are endowed with a self-interaction. In that work the 
spacetime was assumed to be static and spherically sym- 
metric, although the boson stars were only stationary. 

In this paper we want to report that by performing 
a dynamical transition to SC in boson stars where the 
real scalar field 4> is excited to a non-trivial configura- 
tion, we find that the self-interaction term for the boson 
field (other than the mass term) is not required, contrary 
to what was found by Whinnett. This result is not the 
only novelty. In order to perform our analysis we have 
developed several new tools. Unlike Novak's work, we 
always work in the Jordan frame where the interpreta- 
tion of results has a more direct physical meaning. To 
achieve this goal we have constructed a new code based 
on a 3+1 approach for STT developed in [l9[ and adapted 
to a BSSN formulation [2(| [2l|, in which we have imple- 
mented constraint preserving boundary conditions. In 
this regard, it is important to mention that STT have 
been proved to possess a well-posed Cauchy problem in 
the Jordan frame 0, H3- In particular, we shall show 
elsewhere (23j |. that the spherically symmetric equations 
used in this paper are strongly hyperbolic and therefore 
the Cauchy problem is well-posed as well for this partic- 
ular case. 

The paper is organized as follows: In Section [TT| we 
introduce the Scalar- Tensor Theories and discuss briefly 
some properties associated with the Jordan frame. In 
Section [TTT] we describe the boson-star model. Section HVl 
describes the numerical results which include some com- 
ments about the scalar-waves emitted, the details of 
which will appear elsewhere [23]. Finally, Section [V] con- 
tains the conclusions. 



II. SCALAR-TENSOR THEORIES OF GRAVITY 

The general action for STT with a single scalar field is 
given by 

+ S matt [g a b,ip} , (2.1) 

with <f) the non-minimally coupled scalar field, and where 
ip represents collectively the ordinary matter fields (i.e. 
fields other than <fi; we use units such that c = I = Gq). 
For the problem at hand ip will represent the complex 
scalar field that we use to model the boson star (see 

secnni. 

The representation of the STT given by Eq. (I2.1[) is 



called the Jordan frame representation. The field equa- 
tions obtained from the action (|2.1|) are given by 



Gab = 8nT a b 



D(j) 



\f'R = V , 



(2.2) 
(2.3) 



where a prime indicates 8$, □ := <? af, V Q Vb is the covari- 
ant d'Alambertian operator, G a b — Rab — ^gabR, and 
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(2.4) 
(2.5) 
(2.6) 

(2.7) 



where T^ att is the stress-energy tensor for all matter 
fields other that (p. 

Using Eq. (|2.2[) , the Ricci scalar can be expressed in 
terms of the energy- momentum tensor Eq. (|2.4p . Equa- 
tion (12.31) then takes the final form 



□c 



fV - 2f'V - \f (1 + 3p (V^) 2 + \f'T u 
/ (1 + 3/' 2 /2/) 



(2.8) 

where T matt stands for the trace of T^ att . It is important 
to stress the fact that, although we have included here a 
self-interaction potential V(4>) for the STT, for the actual 
analysis of SC in boson stars that we shall present below 
it was not taken into account. 

Using a modified harmonic gauge [19], it has been 
proved that the previous equations (in vacuum) can be 
put in a quasi-linear diagonal hyperbolic form. Moreover, 
such equations can also be put in a full first order 3+1 
form, from which one can prove directly their strong hy- 
perbolic character (22j. Therefore, the Cauchy problem 
is well posed in the Jordan frame. 

The Bianchi identities directly imply 



V r T c 







(2.9) 



However, the use of the field equations leads to the con- 
servation of the energy-momentum tensor of the matter 
alone 



V c T5L t = 



(2.10) 



which implies the fulfillment of the (weak) equivalence 
principle (i.e. test particles follow geodesies of the metric 

g ab -) 



III. BOSON STARS 

As we mentioned in the introduction, the phenomenon 
of SC was first discovered in static neutron star models. 
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Afterwards, the dynamic transition from the unscalarized 
to the scalarized state was also analyzed (l^ . In this 
paper, we propose to study a similar transition but using 
boson stars instead of neutron stars. In some way this 
matter model is even more fundamental than neutron 
stars in that one does not have to assume a perfect fluid, 
but rather a field theory. It is also simpler since one does 
not have to deal with the EOS for the nuclear matter 
(where uncertainties in the model are always an issue), 
as well as the technical difficulties associated with the 
numerical simulation of shock fronts. 

In a previous investigation, Whinnett |17j | constructed 
stationary configurations of boson stars within a STT, 
and showed that SC was only possible if one included 
a self-interaction potential for the boson field. He car- 
ried out the analysis for three different classes of STT, 
one of which in fact is almost identical to the one we 
use in this paper. In our case, instead of constructing 
stationary configurations of scalarized boson stars, we 
analyse the dynamical transition from the unscalarized 
state to the one with scalarization. The final state of this 
process correspond in principle to one of the stationary 
states that one would found using the method by Whin- 
nett. Nevertheless, unlike Whinnett's results, we find SC 
without the need of a self-interaction. The possible ex- 
planation for this will be elucidated in Sec. El We point 
out that boson-star models have been constructed in the 
past within STT [2^ - [27| . but apart from Whinnett's work 
none other study has looked at the spontaneous scalar- 
ization of such objects. 

We shall consider boson stars described by the follow- 
ing stress-energy tensor 



rTimatt 
L ab 



1 " 

2 . 
g a b 



(V„^*)(V 6 V0 + W)(V„^) 



-|VVf + Wl 2 ) 



(3.1) 



where ip is a complex-valued scalar field that rep- 
resents the bosons, |V0| 2 := g ab (V a i/j)(VbilJ*) and 
\ij}\ 2 := tj)*ip. This model arises from the Lagrangian 
Anatt = -\/=9[i|VVf + ^(|Vf)]- The potential 
V^^tp) will be taken to be of the form 



(3.2) 



which includes a mass term but no self-interaction poten- 
tial (which is typically associated with a term hX(ip* ip) 2 ) 

The boson field obeys the Klein-Gordon equation: 



□■0 — m^ip — 



(3.3) 



Notice that the problem we wish to study involves two 
different scalar-fields, a real-valued field <fi coupled non- 
minimally to the curvature, and a complex-valued field 
if) (the boson field) coupled only minimally. 

The matter Lagrangian is invariant with respect 
to a global phase transformation %p —> e iq ip (with 



q a real constant). Noether's theorem then im- 
plies the local conservation of the boson number 
S7 a J a — 0, where the number-density current is given by 
ja = *_ g ab [^v&V** - V*V&# This means that the total 
boson number 



i a J a Vhd 3 x 



(3.4) 



is conserved, where n a is the normal vector to the spatial 
Cauchy hypersurfaces S t and h is the determinant of the 
induced metric on E t . The total boson mass is then given 
by Mbos = rnbAf, with m D '■= hm^/c the mass of a single 
boson (where we have restored the speed of light c and 
the Planck's constant K). 

In fact, when the scalar field (f> acquires a non-trivial 
value one can define another global quantity associated 
to it, which for asymptotically flat spacetimes is given by 



Qscai = lim 



1 

4tt 



(3.5) 



where s a is the unit outward normal to a topological 2- 
sphere S embedded in S t , and r is a radial coordinate 
that provides the area of S asymptotically. 



IV. NUMERICAL RESULTS 

For the analysis at hand, we shall consider a spherically 
symmetric (SS) spacetime with a metric given by 



ds 2 



N 2 dt 2 + A 2 dr 2 + r 2 B 2 (d8 2 + sin 2 9 dip 2 ) , (4.1) 



where the metric coefficients (N, A, B) are all functions 
of the coordinates t and r. The scalar-field variables will 
be functions of t and r as well. As can be seen from the 
form of (|4.ip . we consider a null shift vector. Also, the 
area of 2-spheres is given by A = 4irr 2 B 2 , which only 
coincides with the area coordinates when B = 1. 

We have constructed a spherically symmetric evolution 
code based on the BSSN system of equations, together 
with a 3+1 formulation of the STT equations developed 
in Ref. [l9j]. The details of this system, includin g it s 
hyperbolicity properties, will be reported elsewhere [23j j . 

For the evolution we have also used a generalization of 
the Bona-Masso slicing condition [28| that led to a well 
behaved hyperbolic system [22| . 

We have constructed initial data for stationary boson 
stars like in GR, by assuming ip(t,r) = ^>(r)e lut and 
<f>(t, r) = 0, and solving the eigenvalue problem resulting 
from Eq. (|3.3p using a shooting method to find u. Notice 
that <f>(t,r) = 0, solves exactly Eq. (|2.8[) . The result- 
ing configuration corresponds to a strictly static space- 
time. In order to study the dynamic transition to SC, we 
then considered a small Gaussian perturbation for the 
scalar field 4>. We then solved the new Hamiltonian con- 
straint for this perturbed initial data. The momentum 
constraints are trivially satisfied initially by assuming a 
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moment of time symmetry for which the extrinsic cur- 
vature vanishes, as well as the momenta associated with 
the real scalar field <f> and the boson field ip. 

For the dynamical simulation presented below, we have 
taken the non-minimal coupling function to be of the 
form F(4>) = 87r/(0) = 1 + 87r£0 2 , with £ a positive con- 
stant. 

Several scenarios can happen depending on the initial 
configuration. Figured (top panel) depicts a curve which 
represents the boundary (i.e. the critical values) of the 
transition to SC for different values of the constant £ and 
the central value of the norm of the complex scalar field 
^(0). Initial configurations below the critical line are 
stable with respect to the perturbations and do not lead 
to a SC transition. This means that for such configura- 
tions the scalar field (f> simply radiates away during the 
evolution and leaves behind a stable stationary configu- 
ration with a globally null 4>. On the other hand, initial 
configurations above the critical line (but below the line 
marked "maximum mass") are unstable, and when per- 
turbed lead to a dynamical transition to SC where the 
final configuration is a stationary boson star endowed 
with a non-trivial scalar field <j). 

The horizontal line marked "maximum mass" in Fig- 
ure [T] denotes the maximum mass for stable boson star 
configurations. For boson stars with no perturbation 
from the scalar field <p the mass increases with increas- 
ing \l/(0) until a threshold value vE , (0) cr ;t is reached, af- 
ter which the masses start to decrease for larger ^(0). 
This value separates two regions, the values with ^(0) < 
^(0)crit represent stable boson stars, while the configu- 
rations with ^(0) > x P(0) cr it are unstable and either col- 
lapse to a black hole or migrate to a configuration on the 



stable branch 



The maximum mass is 0.633 and 



corresponds to M r (0) cr i t = 0.076 or er(0) cr it = 0.27 with 
the usual normalization (a = \/47r ty). Similar behavior 
is expected to occur in STT 



Figure [5] depicts a series of snapshots taken at differ- 
ent times during the evolution for the unstable case with 
£ = 6 and ^(0) = 0.03. Note from this figure that ini- 
tially the scalar field (j> almost vanishes, whereas at the 
end of the evolution it settles down into a stationary con- 
figuration which results in a non-trivial profile that inter- 
polates between a finite value (j)(tfi na i,r — 0) at the cen- 
ter of the boson star and a vanishing value 4>{t final i r 00) 
asymptotically. The final stationary configuration is also 
characterized by the appearance of a global scalar charge 
(see Figure Q] bottom panel) . 



V. DISCUSSION 

Boson stars are stable self-gravitating configurations of 
a complex scalar field, and as such they could in prin- 
ciple exist in nature. If the bosons are light (m& := 
hm^/c ~ eV/c 2 ) one expects the mass of the boson star 
to be of order M ~ 10 20 kg, which is in the mass range 
of some asteroids. On the other hand, for heavy bosons 
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FIG. 1: (Top) Critical values of central boson field ^central 
and non-minimal coupling parameter £ for spontaneous 
scalarization; (Bottom) Scalar charge at the end of the tran- 
sition to spontaneous scalarization for different values of the 
parameters. 




FIG. 2: NMC scalar field <j> and metric components at differ- 
ent times for the simulation with £ = 6 and ^(0) = 0.03. 



(nib ~ 100GeV/c 2 ) the boson star mass turns out to be 
much lower (M ~ 10 9 kg). However, by including self- 
interactions it is possible to increase the mass of the bo- 
son star to the order of the solar mass [3(| • It has also 
been speculated that super-massive boson stars instead 
of black holes could be at the center of galaxies [3l| . If 
boson stars actually exist in nature, they could serve as 
natural "laboratories" to test different alternative theo- 
ries of gravity. 

In this work we have used a scalar tensor theory of 
gravity to study dynamical simulations of boson stars. 
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We have found that, just like in the case of neutron stars, 
boson stars can also undergo a spontaneous scalarization 
process. We have analyzed this transition in a dynamical 
fashion using a fully relativistic code in spherical sym- 
metry. Unlike previous studies of stationary boson star 
configurations, we have found that self-interactions are in 
fact not required to produce scalarization. We have im- 
plemented a STT which is very similar to the one studied 
by Whinnett, where scalarization was not found without 
self-interaction. However, he used a value £w = 1 [HI 
(corresponding to our £ = 2) which are actually in the 
region of Figure [T] where scalarization does not ensue. 
He also used £w = 2 (ours £ = 4) where our results of 
Figure Q] show that there is a small region where scalar- 
ization is found. Actually from Whinnett 's figure 4 is not 
clear that scalarization is not found for = 2 since his 
coupling parameter is not null. By increasing the value 
of the non-minimal coupling parameter £ one can reach 
a threshold where the energetically favorable configura- 
tions are those with the presence of a non-trivial scalar 
field 4>. Note from Figure [1] that the larger the value 
of the parameter £, the lower the central energy density 
required to produce the transition to scalarization. 

At this point is perhaps appropriate to men- 
tion that the effective Brans-Dicke parameter given 
by w BD = //(/') 2 Uo (where (f> is the asymp- 
totic (cosmological) value) takes the explicit form 
wbd = (l + 87r^)/[327r(£(/> ) 2 ]. Therefore, as the 
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